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Abstract 

The neutrino emissivity of compact stars is investigated in this work. We 
consider stars consisting of nuclear as well as quark matter for this purpose. 
Different models are used to calculate the composition of nuclear and quark 
matter and the neutrino emissivity. Depending on the model under consider- 
ation, the neutrino emissivity of nuclear as well as quark matter varies over a 
wide range. We find that for nuclear matter, the direct URCA processes are 
allowed for most of the relativistic models without and with strange baryons, 
whereas for the nonrelativistic models this shows a strong dependence on 
the type of nuclear interaction employed. When the direct URCA processes 
are allowed, the neutrino emissivity of hadronic matter is larger than that 
of the quark matter by several orders of magnitude. We also find that the 
neutrino emissivity departs from T 6 behavior when the temperature is larger 
than the difference in the Fermi momenta of the particles, participating in 
the neutrino-producing reactions. 
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I. INTRODUCTION 



The neutrino emissivity of compact stars, such as neutron stars, has been a subject of 
interest for quite some time. When collapsed stellar core is formed after the supernova 
explosion, it initially cools by emitting radiation till the temperature falls to about 10 9 K or 
~0.1 MeV. Below this temperature, cooling by radiation emission is not efficient and the 
star is expected to cool by emitting neutrinos. The neutrino emission occurs by the so- 
called URCA processes ( n — > p + e~ + T7 e and p + e~ — > n + v e ) and the energy-momentum 
conservation requires that the proton Fermi momentum ( p^ ) should be larger than 0.5 xp^. 
This implies that the proton fraction in the neutron star must be larger than 1/9 0. Initially, 
it was thought that such high proton fraction will not be present in neutron stars and one 
has to find some other mechanism of neutrino emission. One of the ways of getting around 
this problem is to consider the so-called modified URCA processes @. In these processes, 
the neutrino emission occurs in the presence of another particle ( a nucleon ) which helps 
in satisfying the energy-momentum conservation. However, the neutrino emissivity due to 
the modified URCA processes is extremely small ||- to explain the cooling rates of these 
stars. Later, it was observed that, if the star contains pion || or kaon condensate J7| or 
if the star consists of quark matter ||- [12], the direct URCA processes are allowed and 
the neutrino emissivity is substantially larger ( ~ order of magnitude ) than that due to 
modified URCA processes. Recently, some authors have found that, for certain models of 
the nuclear matter 13] , the proton fraction could be larger than the critical value above 
which the direct URCA processes are allowed. In such a situation, the neutrino emissivity 
of the star could be large and one can explain the cooling rates in terms of the standard 
neutron star models. This explanation, however, depends on the assumption of the nuclear 
interactions. In particular, it depends on the nature of the three-body force and its isospin 
dependence |L4| . 



The purpose of the present work is to calculate the neutrino emissivity of nuclear and 
quark matter using different models. For nuclear matter, we use nonlinear Walecka model, 
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derivative scalar coupling model, chiral a model and nonrelativistic models, whereas, for 
the quark matter, we use the MIT bag model and the chiral color dielectric (CCD) model. 
One of the reasons behind this calculation is to determine the nuclear equations of state 
in which the direct URCA processes are allowed. This is of some importance, since the 
neutrino emissivity due to direct URCA processes dominates when these are allowed. We 
also want to explore the dependence of the neutrino emissivity on different quark models. 



In our earlier calculations of neutrino emissivity of the quark matter ||17|| , we found that, for 



certain values of Fermi momenta, the approximate neutrino emissivity formula of Iwamoto 



H fails. In that work [17[], we were able to determine an empirical formula which is able 
to reproduce the calculated neutrino emissivity over a wide range of quark matter densities 
and temperatures. In the present calculation, we want to investigate, if a similar situation 
exists for nuclear matter also. 

The results of our calculation can be summarised as follows. We find that for all the 
relativistic models considered here, the direct URCA processes are forbidden if the nuclear 
density is below a certain value. In other words, for all the models, the proton fraction is not 
high enough at lower nuclear densities. The neutrino emissivity, for all the models without 
strange baryons, falls rapidly at higher densities (more than 5 times nuclear matter densities 
). But when the strange baryons are included, the neutrino emissivity due to direct URCA 
process at higher densities varies very slowly with baryon densities. Thus, the appearance of 
strange baryons depletes the neutron fraction and the direct URCA processes are allowed, 
even though the proton fraction is not large. We find that the reactions involving strange 
baryons give a significant contribution to the neutrino emissivity. As for the quark matter, 
we find a large dependence of the neutrino emissivity on the models used, to calculate quark 
matter equation of state. However, when the direct URCA processes in the nuclear matter 
are allowed, the neutrino emissivity of the nuclear matter is larger than that of the quark 
matter at corresponding nuclear density. We also find that the calculated neutrino emissivity 
departs from T 6 [|5|,|T3] behaviour for a certain range of Fermi momenta of the constituents. 
This happens for nuclear as well as for quark matter. Furthermore, we are able to fit the 



calculated neutrino emissivity with a simple universal formula. This clearly implies that the 
departure from T 6 behavior of the neutrino emissivity is kinematical in origin. 

The paper is organised as follows. In Section II, we briefly describe the models used in the 
calculation of nuclear and quark equations of state. In Section III, the neutrino emissivity 
formulae are presented. Finally, the results are discussed in Section IV. 



II. THE MODELS 



A. Nuclear Models 



Four different models, the nonlinear Walecka model ||18|J19| , derivative scalar coupling 



model [2C , 25] , chiral o model [21] and nonrelativistic model [14| have been used to calculate 
the equation of state of the neutron matter. Of these, the first two have been extensively 
used in nuclear structure calculations f22| , |2"0[| and have been able to reproduce the properties 
of nuclei over a wide range of the periodic table. This probably ensures that one has correct 
nuclear equation of state near the nuclear matter density. The chiral a model has been used 
as another model for the nuclear equation of state. This model has been used to calculate 
neutron star properties [51]. One important fact about this model is that nonlinear terms 
can give rise to the three body forces, which is important in the equation of state at high 
densities. For all these relativistic models, the nuclear equation of state is calculated by 
adopting the mean field ansatz. This is in contrast to the nonrelativistic models, where 
actual interaction between the constituents are considered. 

The Nonlinear Walecka Model: (NW) The Lagrangian density of the nonlinear 



Walecka model [|19[] is given by, 
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The Lagrangian in eq.(p])above includes nucleons, A and £~ hyperons (denoted by subscript 
i), electrons and muons (denoted by subscript I) and cr, uj and p mesons (given by a, uj^ and 
pa,n reS p ec tively). The Lagrangian includes cubic and quartic self- interactions of the a field. 
The meson fields interact with baryons through linear coupling and the coupling constants 
are different for nonstrange and strange baryons. 

In the presence of baryons, the mesons develop nonzero vacuum expectation values (a, 
uj and p a respectively). Assuming that the baryon densities are uniform, one finds that the 
time components of uj and p 3 , in addition to a, are nonzero. One can then define effective 
masses (mj) and chemical potentials (/2j) for the baryons as, 

m i = m i - g ai a (2) 

and 

P>i = im- g<j& - hg P NP 3 , (3) 

where J3 is the value of the z-component of the isospin of baryon i. The Fermi momenta 
(ki) and number densities (rii) of the baryons are given by ki = \J pj — mf and rii = 



For leptons, the Fermi momenta and number densities are given by ki = y pj — mf and 
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The parameters of the NW model are meson-baryon coupling constants, meson masses 
and the coefficients of the cubic and quartic self-interactions of a meson ( b and c respec- 
tively). The uj and p meson masses have been chosen to be their physical masses. Of 
the rest of the parameters, the nucleon-meson coupling constants ( — , — and — respec- 
tively ) and the coefficients of cubic and quartic terms of the a meson self interaction ( b 
and c respectively) are determined by fitting the nuclear matter properties ( the binding 
energy/nucleon (— lQMeV ) and baryon density (0.15/m~ 3 ), symmetry energy coefficient 
(32.5MeV ), Landau mass (0.83mAr) and nuclear incompressibility (250 — 300Me\^) ). The 
coupling constants of the hyperon-meson interactions and are not well known. These cannot 
be determined from nuclear matter properties, since the nuclear matter does not contain 



hyperons. Furthermore, properties of hypernuclei do not fix these parameters in a unique 
way. In the literature, a number of choices have been made. These are (a) same as the 



nucleon-meson coupling constants (Universal coupling) pi|,|25|| ,(b) y 2/3 times the nucleon- 



meson coupling constants p^ , p5| , (c) 1/3 times the nucleon-meson coupling constants [18| 
and (d) 2/3, 2/3 and 1 times the nucleon- meson coupling for a, u and p mesons respectively 

We have used the above mentioned choices for meson- strange baryon couplings to in- 
vestigate the neutrino emissivity. 

Derivative Scalar Coupling Model (DSC model): In this case, the Lagrangian is 
given by f20f 

£( x ) = S ^7 M <9 M - nii + ^ayy" - g pi p a ^T a )iJi - ^u^u^ 

i 

+\mlu^ + \{d^a - mla 2 ) - Ifotf + imjpgpg 

+ Y, 9c^Mi{\ + gaCr/M) + - m,)V»i (4) 

i I 

where represents the vector meson field, m a and m w are the masses of the scalar and 
vector fields, and F pv = d p oj v — djAV The summation indices i and I stand for fermions 



and leptons respectively. The DSC Lagrangian differs from the original a - lu model [26 



in the baryon- a meson coupling term (= g^ipipa for the a - uo model), but it has also two 



parameters: g a and g u like the original model [26 



From the Lagrangian (|4]), the following definition of (density dependent) effective nucleon 
mass (m*) is suggested : 

m* = M/(l + g a a/M) (5) 

Chiral Sigma Model (CS model): This model includes a, u and 7r fields. In 
addition, p meson is included in the Lagrangian to reproduce the symmetry energy of the 
nuclear matter correctly. The Lagrangian for an SU(2) x SU(2) chiral sigma model that 
includes (dynamically) an isoscalar vector field (a> M ) is 
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-~£pV< pV^V)V> + ^MWdn ~ m i)ipi ( 6 ) 
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where F^ = c^o^ — d v uj^,%l) is the nucleon isospin doublet, ~7f is the pseudoscalar pion 
field and a is the scalar field. The vector field uj^ couples to the conserved baryonic current 
jfj, — fillip- The expectation value < j > is identifiable as the nucleon number density. 

The interactions of the scalar and the pseudoscalar mesons with the vector boson gen- 
erates a mass for the latter spontaneously by the Higgs mechanism. The masses for the 
nucleon, the scalar meson and the vector meson are respectively given by 

M = g a x 
m a = V2Xx Q 

mu = g^x . (7) 

where x$ is the vacuum expectation value of sigma field. 

With x = (< a 2 + Ix 2 >) 1 ^ 2 , the the effective mass of the nucleon is M* = yM, where 
y = x/xq. For both DSC and CS model, the chemical potential of the baryon can be 
related to the fermi momentum through a similar relation as eq.(|3]). The nucleon- meson 
coupling constant are determined by fitting the nuclear matter properties, binding energy 
(-16.3 MeV), saturation density (0.153 fm~ 3 ) and symmetry energy coefficient (32.5 MeV). 
The incompressibility in the above two models are 225 MeV and 700 MeV respectively. 

Nonrelativistic models: Here we have used the models as proposed by Wiringa et al. 
| 14| by combining different two- nucleon and three- nucleon potentials. In particular three 



different choices have been considered. 

• Argonne v 14 (AV14) and Urbana VII (UVII) three nucleon potential, 

• Urbana w 14 (UV14) two nucleon potential and UVII, 



UV14 and three nucleon interaction (TNI) model of Lagaris and Pandharipande |Tf| . 



AV14 and UV14 have identical structure and can be written as a sum of 14 operator 
components. The main difference between these two models come from the strength of short 
range tensor force. The UV14 does not have short range tensor components which results in 
a weak tensor force that vanishes at the origin. The AV14 tensor force is finite at the origin 
and at intermediate distance looks like Paris potential |L6 |. The three nucleon potential 



potential UVII combines a long range two pion exchange part and an intermediate range 
repulsive part. Out of the three different combinations, AV1A + UVII, UV1A + UVII 
and UV14 +TNI considered, only the second combination produce enough proton fraction 
to have non- zero neutrino emissivity. 

In order to have some idea on the expected behaviour of different models, let us study 
the proton fractions for different parameters as well as different models. In fig. 1, we have 
plotted the proton fraction for nonlinear Walecka model with different hyperon couplings 
and incompressibility. It is evident from this figure that, for all the parameter sets, the 
required proton fraction is attained over a certain density (~ 2p , p is the nuclear matter 
density ). A change in incompressibility from 300 MeV to 350 MeV does not give much 
change in the proton fraction. On the other hand, the variation of hyperon couplings give a 
large change in the proton fraction. This indicates that the effect of hyperon couplings on 
neutrino emissivity will be much larger than that of the incompressibilities. 

Fig. 2. gives a comparison of proton fraction between different models. It shows a strong 
density dependence of the proton fraction on models. Below a certain density (~ 0.2/m~ 3 ), 
none of these models satisfy the limit and neutrino emissivity will be zero. We also find that 
nonlinear Walecka model without hyperons yield larger proton fraction at higher density 
(> 0.8/m~ 3 ). Also, with hyperons required proton fraction is attained at earlier baryon 
densities than without hyperons. However, the nonrelativistic model yield lowest proton 
fraction throughout the density range. 

The proton fraction also depends on the symmetry energy coefficients. We have varied 
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the symmetry energy coefficient from 28- 38 MeV for DSC and CS models. With increase 
in symmetry energy, the required proton fraction for neutrino emissivity in direct URCA 
process, is attained at lower baryon density. 



B. Quark Models 



We have used CCD model and bag model |17j to study the neutrino emissivity. The 



colour dielectric model is based on the idea of Nielson and Patkos ||28|| . In this model, one 
generates the confinement of quarks and gluons dynamically through the interaction of these 
fields with scalar field. Here, we have used chiral extension of this model to study the quark 
matter neutrino emissivity. The CCD model has already been used to study static properties 
of baryons [p9|| , properties of quark matter at finite density and temperature |]30| , |31| , hadron- 



quark phase transition and properties of dense stars ||19| , |32]| . In CCD quark matter, we 
assume that meson (0) expectation value is zero i.e. < <p >— where as < <p 2 >^ 0. So, the 
Lagrangian is rewritten in terms of < <fi 2 >. The quark mass becomes density dependent. 
In this model, with increase in density quark mass decreases and drops to about l/4th of 
its initial value around density 2fm~ 3 . This density is defined as the critical density for 
chiral transition in our model. There are five parameters in CCD model; bag parameter 
B, scalar field potential parameter a, u and d quark mass, strange quark mass and strong 
coupling constant a s . These input parameters are obtained by fitting the baryonic masses. 
In the present paper, we have discussed the results for the parameter set: B l / A = 152.1MeV, 
m q(u,d) = 91-GMeV, m g ( s ) = 294.9MeV r , a = 36 and strong coupling constant ot s = 0.08. In 
the bag model, the neutrino emissivity is calculated with m u = = 0, m s = 150MeV and 
a s = 0.08. For both CCD and bag model, we consider interaction upto first oreder in a s . 

III. EMISSIVITY FORMULAE 

In general, the Lagrangian density for URCA processes in the current current interaction 



form is written as 33 



C(x) = ^Ll^xyL^x) + H.C (8) 
V2 



where the weak coupling constant G F = 1.435 x 10 49 erg cm 3 and 1^ and J M are the leptonic 
and hadronic weak currents respectively. 

l^x) = e7„(l - 7 5 )^ e + /Z7„(l - 75)^ + + /i-c. (9) 

M x ) = ~ + h.c. (10) 

where h.c. stands for hermitian conjugate. For quarks A = B = 1 and for baryons the value 
of A and B depends on the specific nature of the particle. Using the above Lagrangian one 
can calculate the neutrino emissivity e, 

e = gj {n i=1 -0^}E v W fi F{p uP2 , Pe ) (11) 

where i = 1, 2, e, v, E v is the energy of the neutrino and g is the degeneracy factor. The 
transition rate is 

= iMmiMg (12) 

where Pj n is the sum of the initial momenta, Pf n is the sum of the final momenta, Ei is the 
energy of the ith particle and |M| 2 is the squared invariant amplitude averaged over initial 
spins and summed over final spins. The symbol F(pi,p 2 ,p e ) for the reactions 1 — > 2 + e~ + i7 e 
and 2 + e~ — > 1 + u e are 

F d (pi,P2,Pe) =ra(pi)(l-n(p2))(l-n(Pe)) (13) 

and 

F r (Pi,P2,Pe) = (1 - n(pi))n(p2)n(p e ) ( 14 ) 

respectively, where 
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A. Hadronic matter 



In hadronic matter, several weak decays may contribute to the neutrino emission as given 



below 0. 



n — > p + e +z/ e ; p + e — > n + v e 

A = cos6 c ; B = (F + D)cosO c 

A — >■ p + e~ + z/ e ; p + e~ — > A + v e 

A = -3sin9 c /VE] B = - ^ 3F ~L D ^ sin9 c 

v6 

S — > n + e~ + v e \ n + e~ — > S + z/ e 
A = -sin0 c ; B = — (F - D)sin6 c 



A + e~ + v e ; A + e 
A = 0; S 



> E + z/ e 
'^Dcos6 c 



(16) 



(17) 



(18) 



(19) 



where F = 0.427 and D = 0.823. The angle 6 C is the Cabbibo angle and cos# c = 0.948. The 
\M\ 2 for the above processes is given by 

|M| 2 = ^[64(A 2 + B 2 ){(p 2 .p e )( Pl .p u ) + (p 2 .p,)(pi.p e )} 



+64AB{2(p 2 .p e )(pi-Pu) - 2(p 2 .p v )(p l .p e )} 
+64(A 2 - B 2 )m x m 2 {p e .p v )\ 



(20) 



Using the above matrix element, one can calculate the neutrino emissivity for both direct 
and reverse processes. 



e (1^2+e-+;7 e ) 



P2PePv 2 Pidp 2 dp e dp v dp 1 

rmin{ \P2+Pe I, \pi+Pv I } r>2 _i_ 2 _ p2 2 , „2 _ p2 

Ci x / dP (1 - Pl ^ ) (1 + P2 f,* ) 

^maa;{|p2-Pe|,bi— Pv|} 2E\E V 2E 2 E e 



2EiE v 

rmin{\pi+p e \,\p2+Pv\} 2 _i_ 2 _ p2 

+c 2 x / dp (i + ; 

Jmax{|pi-p e |,|p 2 -p^|} 2E 2 E U 



1 



pI+pI-p 2 



—C^mim 2 x 



/ P2 Pe Pv Pidp 2 dp e dp v 
E\E 2 E e 



2E\E e J 
5(E 2 + E e + E v -E 1 )F d ( Pl ,p 2 ,p e ) (21) 
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P2PePu 2 Pidp 2 dp e dp 1/ dp 1 

f rmin{\p2+p e \,\pi+Pv\} r> 2 + n 2 — P 2 Vo + V 2 — P 2 

dx / dP (1 + Pl P » ) (1+ OFF ) 

+C 2 X / dP OF F ) ^ OFF > 

J p 2 PePv Z P\dp2dp e dp v 



—Cj,m\m 2 x 



E\E 2 E e 



S(E 2 + E e -E u -E l )F r (p l ,p 2 ,p e ) (22) 



where Cl = C 2 = and C 3 = 

Similar decays involving instead of e~ are also possible. Here we have considered only 
e~ channel, because e~ density will be much more than the fi~ and hence it will contribute 
more. In URCA processes involving only one hyperon, there is a change in strangeness, so 
the neutrino emission rate which is proportional to sin 2 8 c , is less than one tenth of those for 
nucleonic URCA processes (oc cos 2 6 c ). But in the processes involving only hyperons, there 
is no change in strangeness and hence such reactions are not Cabbibo suppressed. 

The nonrelativistic reduction of the matrix element (eq.(^)) can be obtained by ne- 
glecting the baryon momenta and replacing baryon energy by corresponding masses. The 
reduced matrix element becomes, 

\M\ 2 = ^[64{(A 2 + W 2 )E e E u + (A 2 - B 2 )(p e .p u )} mi m 2 ] (23) 



The corresponding neutrino emissivity formula [[34]] can be derived from eq. (|TT| ) by perform- 



ing the phase space integral using Fermi liquid theory. 

457vr 
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G F 2 C 2 (A 2 + 3B 2 )m 1 m 2 ^ e T b (24) 



B. Quark matter 

In quark matter, the neutrino emission takes place due to following processes. 

d — > u + e~~ + v e \ u + e~ — > d + z/ e , (25) 
s — > u + e~ + u e ; u + e~ — ► d + v e . (26) 
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The matrix element for the quark URCA process is 

|M| 2 d(s) = 6Aco S 2 e c (sin 2 e c )[(p 2 . Pe )( Pl .p u )} (27) 
The neutrino emissivity for the quark URCA processes are given by 

(e ) = C ud (C us ) j PuPePlPd( s )dpudp e dpudpd( s )S(E u + E e - E v - E d ^) 



x 



e (E u -fi u )/T _|_ l e (E e -He)/T _|_ I e {Hd(s)-E d (s))/T _|_ ]_ 

,2 i td2 



^ rmin{\ Pu+Pe \,\ Pd(s)+P „\} ^ ^ ^ p 2 + p 2 _ p2 ^ ^ ^ ^ + p& _ p ^ ^ 
J max{\ Pu - P e\,\ Pd { s )-pA} ^E u E y 2E d ( s )E e 

e d _ >M(s ^ M )(e~) = C du (C su ) J Pd(s)PuPePu 2 dp d ( s )dp u dp e dp u 5(E d ( s) - E u - E e - E u ) 



x 



e (E d(3) -fi d(s) )/T _|_ ]_ e (vu-E u )/T _|_ X e (ne-E e )/T _|_ 1 
x /•mm{b d(s )+P.|,K+Pe|}^ p2 w+ ^_p2 p2 _ p 2 + p 2 

Jmax{\p d ( 3 )-p u \,\p u -pe\} 2E d ( s )E u 2E u E e 

where C ud = C du = 6G ^°f ec , C sn = C MS = ^f^ 1 , G F is the weak decay constant 
and 6*c is the Cabibbo angle. Here and /ij are energy and chemical potential of the i-th 
species (i = u,d,s,e). The single particle energy momentum relation, required to evaluate 
neutrino emissivity for both hadronic and quark matter, is defined in analogy with the 
relation between chemical potential and fermi momentum. 



IV. RESULTS AND DISCUSSION 

In the present neutrino emissivity calculation, we have considered only the direct URCA 
processes for hadron and quark matter. Also, the neutrino emissivity due to direct de- 
cay process and its reverse, at chemical equilibrium, are taken to be equal. The neutrino 
emission from charge neutral hadronic or quark matter at chemical equilibrium, for small 
temperatures (T < lMeV ), occur due to those fermions whose momenta lie close to their 
Fermi surfaces. Therefore, the required kinematical criterion for reaction to occur is that 
the momentum conservation condition be satisfied for the fermion momenta around the re- 
spective Fermi surfaces. For the process 1 — > 2 + e~ + u, the momentum conservation 
condition can be written as p</ + p e F — Pi F = > 0. 

13 



In fig. 3, we have plotted the density dependence of neutrino emissivity at T = 0.5MeV 
from different hadronic models which are nonlinear Walecka model, derivative scalar coupling 
model and chiral sigma model. The quark matter neutrino emissivity using CCD and bag 
model are also plotted in fig. 3. In CCD model, quarks are massive (m( U) d) = !25MeV, 
m s = 300MeV). The u and d quarks masses decrease with increase in density due to non 
zero < 7T 2 >. On the other hand, m s remains unchanged as < K 2 > and < rf > remain 
zero in the medium throughout the range of densities considered |31| . We find that the 



the neutrino emissivity in CCD model is higher (~ order of magnitude) than the the bag 
model, due to the additive effect of quark masses on the neutrino emissivity, the strong 
interaction coupling constant being same (a s is equal to 0.08 for both CCD and bag model 
where a s = g s 2 /4ir). Fig.3. shows that upto density 1.0/m -3 , neutrino emissivity from 
DSC and CS model matter is higher (~ factor of 2- 4) than the other hadronic models. But 
beyond that only nonlinear Walecka model with hyperons dominate. Here one thing should 
be noted that with increase in the baryon density, for all the models, the mean field value of 
Pq increases due to increase in proton fractions. As a result, the neutron energy and hence 
the neutrino energy decreases. After a certain density, neutrino energy becomes negative and 
the reaction stops. This phenomena is more pronounced in case of models without hyperons 
at densities greater than lfm~ 3 . The neutrino emissivity from nonrelativistic model is of 
the same order as the relativistic models, whereas, the neutrino emissivity due to quark 
URCA processes are in general lower compared to nucleonic URCA processes. 

In fig. 4, we have plotted neutrino emissivity from the reaction n — > p + e + u e , with 
temperatures, for two different Ap. The curves for both exact and approximate results are 
given. We find that exact In our calculation we find that our exact neutrino emissivity is 
consistently smaller than the approximate results. In fact, when Ap is large compared to the 
temperature, our numerical result is almost same as the neutrino emissivity obtained using 
the analytic formula (eq. (p4|)). On the other hand, for small Ap(~ T), there is a deviation 
from the analytic result. The similar results were obtained for quark matter system [ITT . 



The deviation of approximate from exact result can be explained in the following ways. 
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In calculating the approximate formula, the neutrino momentum is neglected in the delta 
function. As long as Ap is much larger than T, this approximation gives correct result 
and neutrino emissivity varies as T 6 . On the other hand, when Ap is small (~ T), an 



additional power of T may come in the denominator |T7| and neutrino emissivity is nolonger 
proportional to T 6 . The other reason for the difference comes from the factorization of angle 
and momentum integrals. The momenta may differ from the corresponding Fermi momenta 
by T in the integral. When Ap ~ T, there are reasons in momentum space where Cos9 pn 
(8 pn is the angle between proton and neutron ) is greater than 1, and the rest of the integrand 
is not small. Clearly these regions must be excluded from the integration as these values 
of Cosa are unphysical. If one does not put this restriction, which happens when one 
factorizes angle and momentum integrals, the phase space integral will be overestimated. 

In our earlier section, we have discussed that a strangeness changing reaction is propor- 
tional to Sin 2 8 c . Hence, it would have lower neutrino emissivity compared to the reactions 
where there is no change in strangeness. In the present calculation, we find that the neutrino 
emissivity from n to p decay and £ to A decay are larger (~ order of magnitude)than the 
neutrino emissivity from A to p and E to p decay, which are strangeness changing reactions. 
Overall, the neutrino emissivity from the neutron decay is higher than all the other decays. 

We have studied the neutrino emissivity from hadronic matter for different parameter 
sets. It is found that there is only a small variation in emissivity with incompressibility. 
Also, as mentioned earlier, hyperon couplings are varied in our calculation. With decrease 
in hyperon couplings, it becomes energetically favourable to convert nucleons into hyperons 
as hyperons do not feel the predominantly repulsive force. As a result, with decreasing 
couplings more and more hyperons get populated. This implies that with decrease in hyperon 
couplings, the neutrino emissivity due to the hyperon decay increases. 

In our calculation for neutrino emissivity from hadronic matter, we find that the depar- 
ture from the usual approximation |T7|] arises due to large T/ Ap. The similar behaviour 
have been found for quarks also. So it may be possible to fit the numerically calculated e 
{(-exact) with a function of the form e approa ;.(/(a;))~ 1 , where x = T/Ap. The function f{x) 
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should be such that for small values of x it should approach unity. In fig. 5, we have plotted 
tapprox./texact against T / Ap for different decays with different parameters. We have fitted 
the above graph with a function f(x) = 1 + ax + hx 2 + cx 3 , with a = —2.5, b = 100 
and c = 30 as obtained in the case of interacting quarks ||17| . It is evident that with an 
overall multiplicative factor ~ 1.5, same function can describe both hadronic and quark 
emissivities. This factor 1.5 is due to the nonrelativistic approximation, which is used to 
calculate approximate formula of neutrino emissivity in hadronic decays. 

V. CONCLUSION 

In the present paper, we have studied the neutrino emissivity due to hadronic and quark 
URCA processes for different models. We have also considered effect of incompressibility 
and hyperon couplings on neutrino emissivity. It is found that relativistic models, considered 
here, in general, have nonzero neutrino emissivity and these are higher compared to neutrino 
emissivity due to quark URCA processes. But the scenario is not so simple in case of 
nonrelativistic models. In such models, neutrino emissivity is highly sensitive to the nature 
of three-nucleon interactions. In fact in the present study, we find that only UV14 + UVII 
gives nonzero emissivity. For relativistic models, neutrino emissivity is more sensitive to the 
hyperon couplings than the incompressibilities. 

Our calculation shows that as in the case of quark matter, approximate formula is not 
valid for hadronic weak decays when T j Ap is large i.e Ap is smaller compared to T. An 
alternative formula can be used to calculate the neutrino emissivity in such cases for both 
hadronic and quark matter. 

In conclusion, direct Urea processes in hadronic decays provides an alternate scenario 
for rapid cooling, without the necessity of phase transition to quark matter phase inside 
neutron stars. 
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FIGURES 

FIG. 1. Proton fractions for nonlinear Walecka (NW) model, first three curves are for incom- 
pressibility K = 300MeV 

FIG. 2. Proton fractions for (a) DSC model, (b) CS model, (c) NW model without hyper- 
ons, K = 300MeV, (d) NW model with hyperons, K = 300MeV and (e) nonrelativistic model 
UV14+UVII 

FIG. 3. Variation of neutrino emissivity in c.g.s. units (gm/cm 3 /sec.) with density for 
T = 0.5MeV, (a) DSC model, (b) CS model, (c) NW model without hyperons, K = 300MeF, (d) 
NW model with hyperons, K = 350MeV, (e) CCD model, (f) bag model and (g) nonrelativistic 
model UV14+UVII 

FIG. 4. Neutron decay emissivity in c.g.s. units (gm /cm 3 /sec.) for incompressibility (K) 
300Mey; for ns = 0.4/m~ 3 and Ap = 54.53 (a) Our result (b) Approximate result; for 
tib = 1.4/m -3 and Ap = 0.50MeV (c) Our result, (d) Approximate result 

FIG. 5. Neutrino emissivity ratio is plotted against x = T/Deltap(Deltap = Ap). 
The points corresponds to different hadronic decays. The line corresponds to the function 
f(x) = 1 + ax + bx 2 + cx 3 . a = —2.5 ,b = 100 and c = 30 as obtained in ref.[17] for 
quarks. 
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